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We compute the next-to-next-to-leading order (NNLO) QCD correction to the total cross-section 
for the reaction gg — >■ tt+X. Together with the partonic channels we computed previously, the result 
derived in this letter completes the set of NNLO QCD corrections to the total top pair production 
cross-section at hadron colliders. Supplementing the fixed order results with soft-gluon resummation 
with next-to-next-to-leading logarithmic accuracy we estimate that the theoretical uncertainty of 
this observable due to unknown higher order corrections is about 3% at the LHC and 2.2% at the 
Tevatron. We observe a good agreement between the Standard Model predictions and the available 
experimental measurements. The very high theoretical precision of this observable allows a new 
level of scrutiny in parton distribution functions and new physics searches. 
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INTRODUCTION 

Production of top quark pairs at hadron colliders is 
among the processes that are most challenging to the- 
ory. Bringing this process under good theoretical control 
therefore represents a significant step in our broader un- 
derstanding of perturbative Quantum Chromodynamics 
(QCD) and its applications at hadron colliders. 

The first step in this direction was made some 25 years 
ago, when the next-to-leading order (NLO) QCD correc- 
tions to tt production were computed in the groundbreak- 
ing works Q- The complexity of the NLO calculations 
required the application of purely numerical methods, 
and it took almost twenty years before the exact analytic 
result appeared Q revealing the full complexity of the 
cross-section for massive fermion hadroproduction. 

In the last few years we are witnessing a signifi- 
cant interest in computing next-to-next-to leading order 
(NNLO) corrections to hadron collider processes. Such a 
demand is dictated in part by the high-precision mea- 
surements available from the LHC and the Tevatron. 
The first hadron collider processes that were computed at 
NNLO, namely, Drell-Yan and vector boson [4-6], Higgs 
and diphoton [l^] production, all share the proper- 
ties of (a) having massless QCD partons and (b) involv- 
ing at leading order (LO) two partons meeting in a color 
singlet vertex. Tackling processes with higher complex- 
ity, among which tt production is a prominent example, 
proved to require new computational approaches. 

About one year ago, the first step in this direction was 
made precisely in the context of tt production. Based on 



a new view 



ll| about how to treat double-real radiation 



corrections, the first genuinely NNLO corrections to the 
total inclusive cross-section in qq ^ tt + X were com- 
puted Later on, the partonic reactions involving 



at least one fermion in the initial state were also com- 
pleted IJ, llJ] . In this work we report the calculation 



of the last missing NNLO correction to tt production, in 
the partonic reaction gg — > tt + X. With this calcula- 
tion, the complete set of NNLO corrections to the total 
inclusive cross-section for top pair production at hadron 
colliders is now known. In this letter, for the first time, 
we quantify their phenomenological implications. 

Before closing this section we would like to point out 
the very recent NNLO calculation of the process pp 
H + j [l5| which was performed with methods similar to 
ours and, in par ticular, the subtraction scheme proposed 
by one of us Moreover, a first partial result for dijet 
production pp — >• jj at NNLO has just appeared We 
believe that this burst of precision applications at hadron 
colliders marks the outset of a new and lasting stage in 
precision physics at hadron colliders. 



THE tt PRODUCTION CROSS-SECTION 

In this letter we consider the total inclusive tt produc- 
tion cross-section 
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The indices i,j run over all possible initial state par- 
tons; /3max = \/l — 4m2/5; is the cm. energy of the 
hadron collider and /? ~ \f\~~p, with p = Am? / s, is the 
relative velocity of the final state top quarks with pole 
mass m and partonic cm. energy ^/s. 

The function $ in Eq. ([T]) is the partonic flux 
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expressed through the usual partonic luminosity 
£ij{x,np) =x{f,(g}fj){x,fi%). 
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As usual, /K_R,_F are the renormalization and factoriza- 
tion scales. Setting jjLp = = m, the NNLO partonic 
cross-section can be expanded through NNLO as 
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In the above equation as is the MS coupling renormalized 
with Nl — 5 active flavors at scale = and (t|"^ 
are functions only of /?. The procedure for restoring the 
dependence on ^ ^ m is standard and has been 
detailed, for example, in Ref. 14 1. 

All partonic cross-sections are known exactly through 
NLO [iHsl. The NNLO corrections to the partonic reac- 
tions ij = qq, qg, qq, qq' , qq' were computed in Refs. [12- 
14j . In the following we present the results for ij = gg. 
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PARTON LEVEL RESULTS FOR gg^tt + X 



FIG. 1: The functions Fj*"'], (|9I10I11|) versus the 80 computed 

points Fi — Fl^\i = 2,1,0 (numerical errors at each point 

are also shown) . For improved visibility, the function f[^^^ is 

(flt) 



Keeping the dependence on the number of light flavors muhiplied by a factor of 10, while * by a factor of 10*. 
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Nl explicit, the NNLO correction agg reads 



The functions Fi = F, 
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- 0, 1, 2, read: 
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Fi = 4gM(-0-00611924+0.0436508L/3)/;3 
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and with the constant Cgg = (as defined in Ref. .17]). 

The functions ^2,1,0 are computed numerically, in 80 
points on the interval /? G (0,1). Details about the cal- 
culation are given in the next section. 

Following the approach of Refs. 
(fit) 



12|-|l4j. the functions 
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29.1805L^ 

-h(-0.0240072 -I- 1.81537L^ -f- 3.14286L^)//3](9) 
10"'' [(6.44022;3 - 4.8664;32 _ o.0324653L2)p 
-h(-13.8424;3 -I- 4.7366;32 - 2.91398Lp)p2 
+(8.43828/3 - 2.78748/3^-1- 2. 38971;33)p3] , (9) 

-0.0195046/3 - 1.4717/3^ - 0.223616/3^ 
+0.499196/3^ + 1.32756;9^ + 0.00466872/33l^ 
+0.0321469/3^L^ + (0.579781L2 + 0.166646L3)p 
+ (-1.36644Lp + 2.24909L^)p2^ (10) 

581.27542;3 + 1251.4057/3^ - 60.478096/3^ (11) 
+1101.2272/3^ - 2905.3858/3^ + 629.9128/3% 
-5.1891075Lp + (1200.741Lp + 162.50333L2)p 
+ (36.074524Lp - 1192.8918L2 - 1810.2849/3)^^ 
+ 1568.7591/3p^ - 461.21326;9p^ + 121.6379/3^^ , 



F2 I Q are derived as fits to the difference Fi 
The functions F^-^^^ together with the discrete values for 
Fi — Ff"^ (including the numerical errors) are shown in 
fig. [1] As can be seen from fig. [T] the functions Ff^^'^ 
vanish smoothly at threshold ;8 — >■ 0, which implies that 
our calculation agrees with the exactly known threshold 
behavior [l3|. This is a strong check of our result. 

To assess the size of the newly derived NNLO correc- 
tion, in fig. [2] we compare: (a) the exact NNLO result, 
(b) the approximate NNLO result with exact Born term 
and (c) the approximate NNLO result with Born term 
restricted to its leading power of /3. Each of these three 
partonic cross-sections is multiplied by the gg partonic 
flux Eq. dll) for LHC 8 TeV. We observe that the power 
corrections derived in the present work are very large. 
In fact their contribution to the integrated cross-section 
is virtually as large as the one due to pure soft gluon 
corrections. 

The partonic cross-section's leadin g po wer behavior in 
the high-energy limit ^ — > 1 reads [H Il8l - E^ 



where Lp = ln(p) and L^ = ln(/3). The functions K 
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constitute the analytically known threshold approxima- 
tion to a"gg [TtI , including the exact Born term 
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The constant c\ « —5.1891075 ... is known exactly [23| . 
To improve the accuracy of the partonic result ([5]) in the 
high-energy limit, we have imposed on it the logarithmic 
behavior ~ c\ ln(p) implied by Eq. (jl3l) . Numerical pre- 
diction for the constant term cq was given in Ref. 24 1. 
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FIG. 2: Partonic cross-section times gg flux ((2]) for the follow- 
ing three cases: exact NNLO (thick black line), approximate 
NNLO with exact Born term (blue dashed line) and approxi- 
mate NNLO with leading Born term (thin red line). 



Our fits return the value cq = —31.96 -I- 0.1119A''i which 
falls within the range estimated in Ref. 24 1. 

The parton level results derived in this section can be 
used to derive an estimate for the so-far unknown con- 
stant Cg'f appearing in the threshold approximation . 



Expanding Eq. [5] around the limit /3 — > we obtain 



C^f = 338.179- 



26.8912iVi, + 0.142848A^£ . 



(14) 

As explained in Ref. [25|, the estimate ([T4| for Cgg 
has to be used with caution and a sizable uncertainty 
should be assumed. We have no good way of estimating 
the error on the extracted constant and to be reasonably 
conservative in the following we take this error to be 50%. 

(2) I li 

The constant Cgg is related [26[ to the hard matching 

(2) 

coefficients Hgg\ g needed for NNLL soft gluon resumma- 
tion matched to NNLO. However, since our calculation 
deals with the color averaged cross-section, we cannot 

(2) 

extract both constants Hgg\ g. We proceed as follows. 
Close to threshold, the color singlet and color octet 



contributions to a 



(2) 



jgg have independent constant terms 

(2) (2) 

^gg,i,S' '^it^ constant Cgg 



in Eq. ([T4| being their 
color average. We parameterize the second, unknown, 



(2) (2) 

combination of Cgg\ g by their ratio Rgg 
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which has the advantage of being normalization inde- 

(2) 

pendent. For any guessed value of Rgg , together with 
Eq. we can extract values for the hard matching 

As a guide for a reasonable value of 



constants H, 

?(2) 



(2) 

99,1,8- 



R^g we take the one-loop result (see [17|, l25|): R 



Cf^sicfl, = 2.18. 
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In the following we vary Rgg in the range 0.1 < R\^g < 

(2) 

8; for each value of Rgg we then vary the color av- 

(2) 

eraged constant Cgg' by additional 50%. We observe 



?(2) 



that as a result of this rather conservative variation, 
the NNLO-FNNLL theoretical prediction for LHC 8 TeV 
changes by 0.4% (in central value) and by 0.2% (in scale 
dependence). Given the negligible phenomenological im- 
pact of these variations, we choose as our default values: 

^Si = 53-17, ^^8 = 96.34 (for Nl ^ 5) , (15) 

derived from Eq. (|14l) and the mid-range value Rgg = 1. 



CALCULATION OF gg tt 



X THROUGH 0{as) 



The calculation of the 0{ag) corrections to gg 



tt 



X is performed in complete analogy to the calculations 
of the remaining partonic reactions |12|-|14|. The two- 



loop virtual corrections are coniputed in [27| . utilizing 
the analytical form for the poles [28j . We have computed 
the one-loop squared amplitude; it has previously been 
computed in [29| . The real- virtual corrections are derived 
by integrating the one-loop amplitude with a counter- 
term that regulates it in all singular limits 'SO^. The 
finite part of the one-loop amplitude is computed with 
a code used in the calculation of pp tt + jet at NLO 
31 1 . The double real corrections are computed in jll[ . 
Factorization of initial state collinear singularities as well 
^J'F,R scale dependence is computed in a standard way; 
see Refs. ^EM- 



PHENOMENOLOGICAL APPLICATIONS 

In table |T] we present our most precise predictions 
for the Tevatron and LHC at 7, 8 and 14 TeV. 
All numbers are coniputed for m = 173.3 GeV and 
MSTW2008nnlo68cl pdf set [12] with the program 
Top++ (v2.0) fsi]. Scale uncertainty is determined 
through independent restricted variation oi and /i^,. 
Our best predictions are at NNLO and include soft gluon 



Collider 


cTtot [pb] 


scales [pb] 


pdf [pb] 


Tevatron 


7.164 


+0.110(1.5%) 
-0.200(2.8%) 


+0.169(2.4%) 
-0.122(1.7%) 


LHC 7 TeV 


172.0 


+4.4(2.6%) 
-5.8(3.4%) 


+4.7(2.7%) 
-4.8(2.8%) 


LHC 8 TeV 


245.8 


+6.2(2.5%) 
-8.4(3.4%) 


+6.2(2.5%) 
-6.4(2.6%) 


LHC 14 TeV 


953.6 


+22.7(2.4%) 
-33.9(3.6%) 


+ 16.2(1.7%) 
-17.8(1.9%) 



TABLE I: Our best NNLO+NNLL theoretical predictions for 
various colliders and cm. energies. 



resummation at NNLL 
In this letter we take A 



as a default value for the 
The reason for switch- 



35[ 



constant A introduced in Ref. 
ing to a new default value for A (compared to A ~ 2 in 
12 -[13, 26 1) is that this constant is consistently defined 



only through NLO. Nonetheless it contributes at NNLO 



4 



too, and a consistent NNLO treatment would require the 
analysis of Ref. 35| to be extended to NNLO, which is 
now possible with the help of the results derived in this 
letter as well as Ref. Given the numerical effect is 

small (a 0.7% shift at LHC 8 TeV and a 0.4% shift at the 
Tevatron), in this work we take A ^ 0. 

As can be concluded from table |T] the precision of the 
theoretical prediction at full NNLO+NNLL is very high. 
At the Tevatron, the scale uncertainty is as low as 2.2% 
and just slightly larger, about 3%, at the LHC. The inclu- 
sion of the NNLO correction to the gg-initiated reaction 
increases the Tevatron prediction of Ref. [13] by about 
1.4%, which agrees well with what was anticipated in 
that reference. 



Collider 


o-tot [pb] 


scales [pb] 


pdf [pb] 


Tevatron 


7.009 


+0.259(3.7%) 
-0.374(5.3%) 


+0.169(2.4%) 
-0.121(1.7%) 


LHC 7 TeV 


167.0 


+6.7(4.0%) 
-10.7(6.4%) 


+4.6(2.8%) 
-4.7(2.8%) 


LHC 8 TeV 


239.1 


+9.2(3.9%) 
-14.8(6.2%) 


+6.1(2.5%) 
-6.2(2.6%) 


LHC 14 TeV 


933.0 


+31.8(3.4%) 
-51.0(5.5%) 


+16.1(1.7%) 
-17.6(1.9%) 



TABLE H: Pure NNLO theoretical predictions for various 
colliders and cm. energies. 



To assess the numerical impact from soft gluon re- 
summation, in table |TT] we present results analogous to 
the ones in table |I] but without soft gluon resummation, 
i.e. at pure NNLO. Comparing the results in the two 
tables we conclude that the effect of the resummation 
is a (2.2, 2.9, 2.7, 2.2)% increase in central values and 
(2.4, 2.2, 2.1, 1.5)% decrease in scale dependence for, re- 
spectively, (Tevatron, LHC7, LHC8, LHC14). 

Next we compare our predictions with the most precise 
experimental data available from the Tevatron and LHC. 
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FIG. 3: Theoretical prediction for the Tevatron as a function 
of the top quark mass, compared to the latest combination of 
Tevatron measurements. 
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FIG. 4: Theoretical prediction for the LHC as a function of 
the collider cm. energy, compared to available measurement 
from ATLAS and/or CMS at 7 and 8 TeV. 



The comparison with the latest Tevatron combination 
[36| is shown in fig. |31 The measured value txtot = 7.65 ± 
0.42 pb is given, without conversion, at the best top mass 
measurement 



37| 



m = 173.18 ± 0.94 GeV. From this 
comparison we conclude that theory and experiment are 
in good agreement at this very high level of precision. 

In fig. m we show the theoretical prediction for the 
tt total cross-section at the LHC as a function of the 
cm. energy. We compare with the most precise avail- 
able data from ATLAS at 7 TeV 0, CMS at 7 ^ and 
8 TeV as well as the ATLAS and CMS combination 
at 7 TeV [4l|]. We observe a good agreement between 
theory and data. Where conversion is provided [s^ , the 
measurements have been converted to m = 173.3 GeV. 

Finally, we make available simplified fits for the top 
mass dependence of the NNLO+NNLL cross-section, in- 
cluding its scale and pdf uncertainties: 



a{m) = (j{mref ) ( ^^^^^ ) 

V TO / 



1 + ai 



TO 

TO — rriref 



(16) 



TO 



ref 



0.2 



771 — rriref 
rUref , 



The coefficient ai_2 can be found in table Hill 



CONCLUSIONS AND OUTLOOK 

In this work we compute the NNLO corrections to 
gg tt + X . With this last missing reaction included, 
the total inclusive top pair production cross-section at 
hadron colliders is now known exactly through NNLO 
in QCD. We also derive estimates for the two-loop hard 
matching coefficients which allows NNLL soft-gluon re- 
summation matched consistently to NNLO. All results 
are implemented in the program Top++ (v2.0) [33| . 
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ruref = 173.3 GeV 


aijnref) [pb] 


ai 


a2 




Central 


7.1642 


-1.46191 


0.945791 




Scales + 


7.27388 


-1.46574 


0.957037 


Tevatron 


Scales — 


6.96423 


-1.4528 


0.921248 




PDFs + 


7.33358 


-1.4439 


0.930127 




PDFs - 


7.04268 


-1.4702 


0.936027 




Central 


172.025 


-1.24243 


0.890776 




Scales + 


176.474 


-1.24799 


0.903768 


LHC 7 TeV 


Scales — 


166.193 


-1.22516 


0.858273 




PDFs + 


176.732 


-1.22501 


0.861216 




PDFs - 


167.227 


-1.2586 


0.918304 




Central 


245.794 


-1.1125 


0.70778 




Scales + 


252.034 


-1.11826 


0.719951 


LHC 8 TeV 


Scales — 


237.375 


-1.09562 


0.677798 




PDFs + 


251.968 


-1.09584 


0.682769 




PDFs - 


239.441 


-1.12779 


0.731019 



TABLE III: Parameters of the fit ([T6| to our best 
NNLO+NNLL prediction, including scale and pdf uncer- 
tainty, for the top pair cross section at the Tevatron and LHC 
at 7 and 8 TeV. The fits are accurate to within few permil in 
the 130-210 GeV mass range. 



The theoretical precision achieved in this observable is 
very high. To illustrate this, we compare to a NLO level 
calculation. At LHC 8 TeV we observe a decrease in scale 
dependence by a factor of (4.3, 4.2, 3.0) when compared 
to, respectively, (NLO, NLO-I-LL, NLO-I-NLL). The cor- 
responding numbers for the Tevatron are (3.9,4.1,2.0). 

The predicted ti cross-section agrees well with all avail- 
able measurements. We are confident that its very high 
precision will enable a new generation of precision col- 
lider applications to, among others, parton distributions 
and searches for new physics. 
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